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Abstract 

In this article, we present a graphical user interface software for the litter de- 
composition and soil carbon model Yasso07 and an overview of the principles 
and formulae it is based on. The software can be used to test the model and 
use it in simple applications. Yasso07 is applicable to upland soils of different 
ecosystems worldwide, because it has been developed using data covering the 
global climate conditions and representing various ecosystem types. As input 
information, Yasso07 requires data on litter input to soil, climate conditions, 
and land-use change if any. The model predictions are given as probability 
densities representing the uncertainties in the parameter values of the model 
and those in the input data - the user interface calculates these densities 
using a built-in Monte Carlo simulation. 

Keywords: decomposition, Monte Carlo simulation, software, soil carbon, 
statistical modelling, uncertainty estimation 



"YassoO? modelling project: www. environment. fi/syke/yasso 

*The corresponding author. 
Email addresses: mikko.tuomi@utu.fi (M. Tuomi), jussi.rasinmaki@simosol.fi 
(J. Rasinmaki), jari.liski@ymparisto.fi (J. Liski) 

^Finnish Environment Institute, Ecosystem Change Unit, Natural Environment Centre, 
Mechelininkatu 34a, P.O.Box 140, 00251, Helsinki, Finland 

^University of Helsinki, Department of Mathematics and Statistics, Gustaf Hallstromin 
katu 2b, P.O.Box 68, Helsinki, Finland 

■^Simosol Oy, Etelainen Asemakatu 2b, 11130, Riihimaki, Finland; Telephone: +358 44 
040 5859 

^Telephone: +358 40 748 5088; Fax: +358 20 490 2390 



Preprint submitted to Environmental Modelling & Software 



January 25, 2013 



Software availability 

Software name: Yasso07 User Interface (Yasso07-UI) 

Developer: J. Liski, Finnish Environment Institute 

Software engineer: J. Rasinmaki, Simosol Oy 
5 First available: 2009 

Hardware requirements: None 

System requirements: Windows, Linux, or OS-X 

Software requirements: Python 2.5 on OS-X 

Availability: www.environment.fi/ syke/yasso 
10 Source code: code. google. com/p/yasso07ui 

Programming languages: Fortran 90 (Yasso07), Python (Yasso07-UI) 

Software package size: 39Mb (Windows), 34Mb (Linux/OS-X) 



1. Introduction 



The decomposition of organic matter (OM) is, in addition to photosyn- 
thesis, the other important process that regulates the terrestrial carbon cycle. 
Decomposition controls carbon dioxide (CO2) emissions from soils into the 
atmosphere and affects the carbon stocks in the soils. Estimates of these 
emissions and the stock are therefore needed in order to estimate the terres- 
trial carbon budget. Recently, concern about climate change has increased 
interest in terrestrial carbon cycle. Terrestrial ecosystems act as remarkable 
sinks and sources of atmosphe ric CO2 and, consequently, have a significant 
effect on climate fllPCCl . l2007h . 

The carbon cycle of soils has been modelled using a variety of a pproaches 
and techniq ues, and diffe rent soil carbon models exist; e.g. Century flPart on et al. 



1983,ll992|), Cou pModel (IJansson and Karlberd. 120041) . Q-model f Rolff and Agrenl . 



19991 ). RO MUL (IChertov et all. 1200 11) . RothC fIColeman and Jenkinsson 



DEC OMP (Wallman et al.l .l2006h. and Yasso07 fiLiski et al.LboOSHTuomi et al 



20051 



Post et al 



20091. l2010a). Du e to the high spatial variability of soil cabon stocks (e.g . 



19821 ) and high uncertainty in their changes (e.g. iPost et al 



2OOII ). and because measuring these variables directly is difficult, laborous, 



and expensive, soil carb on models are commonly used to estimate these stock s 
and their changes (e.g. iPeltoniemi et al.l . l2006l 120071 : iMakipaa et al.l . |2008[ ). 

Some of the aforementioned models describe soil carbon cycle at a rather 
detailed level. For this reason, they also require detailed input information. 
Such information is not always available, which makes it difficult to apply 
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these models in large geographical scales or at national level. However, the 
Yasso OT model requires only litt le easily accessible input information to op- 
erate fiTuomi et al.l . boool . l2010ah . 



S oftwar e are already available to operate the CENT URY ( jParton et al 



19871 . Il992[ ) and RothC ( IColeman and Jenkinssoru . |2005| ) soil carbon models 



but the growing number of practical applications has increased the demand of 
easily accessible soil carbon models (). Acknowledging this demand, we con- 
structed a graphical user interface software for the YassoO? model. This soft- 
ware, hereafter Yasso07-UI, can be used to operate the model and study its 
performance in estimating the litter decomposition and soil carbon cycle prior 
to possibly implementing the model into other simulation softwares. Yasso07- 

UI is also s uitable for various practical applications, such as greenhouse gas 

inven tories (IFinland's Fifth National Communication under the United Nations Framework Conve 



20091 ) and estimation of the effects th e removal of forest harvest residues has 



50 on carbon balances of boreal forests (iRepo et al.l . l201ll ) 



Estimating the modelling uncer tainties as reliably as possible has been a 



topic in several recent studies (e.g. Post et al. . 2008 : Updegraff et al. . 2010l : 



O'Haganl . l201ll ). We acknowledge the importance of accounting for all the 
sources of uncertainty. In this work these sources are (1) the uncertainties 
of the model parameters; (2) the uncertainties in the values of the input 
variables, such as the amount of OM initially in the soil and the litter input 
to the soil; and (3) the uncertainties in the chemical quality of the litter 
input. Our software takes all these sources of uncertainty into account by 
using a built-in Monte Carlo simulation. 



2. Model structure 

Yasso07 model describes litter decomposition and soil carbo n cycle based 



on th e chemical quality of the OM and climatic conditions ( iTuomi et al 



20091 ). Decompos ition of woody litter depends additionally on the physical 



size of the litter (ITuomi et al.l . l2010al ). The model works by dividing fresh 
OM, e.g. leaf, fine root, and woody litter, into four chemically distinguish- 
able fractions that decompose at their unique rates. These fractions are water 
solubles (W), ethanol solubles (E), acid hydrolysables (A), and compounds 
neither soluble nor hydrolysable (N). In addition, there is a humus (H) frac- 
tion, assumed to consist of more recalcitrant compounds, that receives a part 
of the decomposition products of the A, W, E, and N fractions. Carbon flows 
between the different fractions are shown in Fig. [1] (for detailed mathemat- 



3 



ical f ormulae describing the Yasso07 model structure, see also iTuomi et al.l . 
2009h . 

Wit h one exception, Ya s so07-U I is based on the same structure as was 
used in iTuomi et al.l ( 120091 . l2010a( ). Since there were systematical differ- 
ences in the mass loss rates between the European and American litter bag 
measurements used to determine the pa rameter values of the model earlier 



( jPalosuo et al.l . l2005t iTuomi et al.l . l2009l ) . we modelled these differences in a 



more consistent manner by assuming that a fraction of the OM in the litter 
bags is transferred out of the bags with varying rates. This process called 
leaching, can be driven by water flowing through the litter bags. Since the 
mesh size was different between the European and American litter bags (1.0 
mm and 55 fim, respectively), it is natural that the leaching rates can also 
be different. The different mesh sizes have an effect also by filtering faunal 
decomposers in different manners. We modelled the different leaching rates 
as additional free parameters and replaced the original model (Eq. (1) of 



fiTuomi et al.l . 12009])) by 



X(t) = Ax(t) + b - UJ.lPa, X(0) = Xo, 



(1) 



where x is a vector describing the amount of organic carbon in each mod- 
elled compartment, b is the input vector to the system, A is a decompo- 
sition matrix, Ui,i = 1,2 are free parameters describing the precipitation 
induced leaching rates in European (oJi) and American {U2) litter bags, 
I = (1, 1, 1, 1, 1)"^ is a constant column vector, and Pais the annual precipita- 
tion. We selected this model formulation because this effect was not present 
exclusively in the water soluble compartment W. We tested a model with 
I = (0, 1, 0, 0, 0)^, where the unity corresponds to the compartment contain- 
ing the water soluble compounds, but this model was found to have much 
lower posterior probabihty and it could not fully account for the differences 
between the European and American litter bags. 

The far greater mesh size of the European litter bags likely enables greater 
leaching rate. This difference may also have an effect on the faunal acces- 
sibility to the litter inside the bags, which co uld increase the r nass loss rate 
in lit ter bags with greater mesh size slightly ( iBradford et al.l . l2002l : llrmlerl . 
2000). However, this effect could not be quantified from the available data. 
Therefore, we simply analysed the data and estimated the parameter prob- 
ability density functions (PDF's) using Eq. ([T]) but calculate the model 
predictions in the Yasso07-UI using the original equation because the OM 
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Figure 1: Flow diagram of the Yasso07 model and the relative magnitudes of each mass 
flow between labile compund groups of organic carbon and more calcitrant humus; acid hy- 
drolysables (A), water solubles (W), ethanol solubles (E), and compounds neither soluble 
nor hydrolysable (N). The carbon flows whose magnitudes differ statistically (95% confi- 
dence) from zero (solid arrows) between and out of the A, W, E, and N fractions (square 
boxes); the small flows (dotted arrows) into humus (bottom box), each approximately 
0.5%; and the mass flows (dashed arrows) whose MAP estimates were indistinguishable 
from zero but whose 95% Bayesian confidence interval was broader than 0.05. See Table 
[l]for details. 
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may be leached out of the htter bags but not out of the whole soil system. 
The parameters 0Ji,i = 1,2 can in fact be considered as only nuisance pa- 
rameters that help understanding the process of making the measurements 
but not the process of decomposition. This impr ovement had minor effects 
on the model parameter values - see Table [T] and iTuomi et al.l (120091 ). 

The e xtension of the mod el into the woody litter regime is that intro- 
duced in (ITuomi et al.l . l2010al ). Its three parameters describing the effect the 



physical size of pieces of woody litter have on the decomposition, namely, (f)i, 
02, and r,are also shown in Tabled! 



3. The coverage of measurements and model predictions 

The data used to construct the Yasso07 model consisted of several indi- 
vidual datasets describing the mass lo ss rate of leaf, fine root, and woody 
litter in a variety of climatic conditions ( ITuomi et al.l . l2009l . l2010al ). The time 
and climate coverage of measurements is the range of reliable applicability 
of the Yasso07 model and the Yasso07-UI introduced here. Despite the fact 
that the model could be applicable beyond this coverage, the reliability of 
such applications cannot be known. Therefore, we describe this range briefly. 

The climatic range of the measurements of non-woody litter covers more 
than 90% of the globally availa ble climatic conditi ons on land, excluding only 
the driest deserts and glaciers (ITuomi et al.l . l2009l ). The litter types included 
in the data used to develop the model in terms of the initial chemical compo- 
sitions of litter, cover a wide variety of litters from conifers and broadleaved 
trees to shrubs and grasses, covering also a variety of different initial nitrogen 
concentrations (ITuomi et al.l . l2009l ).The measurements used to develop the 
model include an extensive data set on decomposition of non-woody litter 
across Europe, and North and Central America (A^ = 9605), data sets on the 
decomposition of woody li tter in Finland and n eighboring regions in Estonia 
and Russia {N = 2102) (iTuomi et all . l2010ah . Therefore, the Yasso07-UI 



can be applied reliably to a wide variety of litters from different species in 
almost any climatic conditions. 

The time-span of non- woody litter mass-loss measurements covers the first 
twelve years of decomposition of fresh litter extensively. In longer timescales, 
of hundreds of years, the litter becomes less important and the decomposi- 
tion of humus starts to dominate the model predictions because of its low 
decomposit ion rate. These timescales are also covered by the available mea- 
surements (iLiski et al.l . Il998l . l2005l : iTuomi et al.l . |2009| ) and the model is rea- 
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Table 1: MAP parameter values of Yasso07 and their 95% confidence limits. 

Parameter Value Unit Interpretation 

aA 0.72±0.09 a~^ decomposition rate of A 

aw 5. 9 ±0.8 a~^ decomposition rate of W 

ue 0.28lo;o4 a"-*^ decomposition rate of E 

CKjv 0.0311o;qo4 a"^ decomposition rate of N 

pi 0.48±0.06 - relative mass flow, W A 

P2 O.Ol;l:o;oi - relative mass flow, E A 

P3 0.83^0:23 - relative mass flow, N -> A 

P4 0.99to:o5 - relative mass flow, A ^ W 

P5 0.00+0 00 - relative mass flow, E ^ W 

Pe 0.01+0 0? - relative mass flow, N W 

P7 0.00^0;°,') - relative mass flow, A E 

Ps 0.00+o;oo " relative mass flow, W E 

pg 0.02to;o2 - relative mass flow, N ^ E 

Pio O.OO+oioo - relative mass flow, A ^ N 

pii 0.015±0.015 - relative mass flow, W ^ N 

Pi2 0.95lo:5'g - relative mass flow, E N 

cui -0.151±0.008 a~^m~^ precipitation induced leaching (Europe) 

UJ2 0.000+0 002 a^-'^m"^ precipitation induced leaching (Americas) 

Pi 9.5±2.0 10~^ °C~^ temperature dependence 

^2 -1-4+09 -'-0"^ °C~^ temperature dependence 

7 -1.21±0.14 m~^ precipitation dependence 

Ph 4.5±0.8 10~^ mass flow to humus 

an l-6la2 10~^ a~^ humus decomposition rate 

01 -1.71±0.16 cm~^ first order size dependence 

02 0.86 ±0.10 cm~^ second order size dependence 
r 0.306±0.013 - size dependence power 
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sonably trust worthy in predictin, 



scale as well ( iTuomi et al.l . 12010 



the carbon stocks of soils in the global 
This makes the model trustworthy in 
predicting the mass loss of non-woody litter in all the timescales. 

For woody litter, the measurements used to develop the model include 
branches and stems ranging from 0.5 to 60 cm in diameter, and the mass 
loss of these woody biomass con iponents has been f ollowed for 1-70 years 



since the start of decomposition (ITuomi et al.l . l2010a] ). Hence, the model is 



150 reliable in describing decomposition of wood in this timescale. 



4. YassoOT User Interface 

The Yasso07-UI package contains the YassoOT model core and the YassoOT- 
UI software that exploits this model when calculating the predictions. There 
is also a sample of 100 000 vectors drawn randomly from the posterior prob- 

155 ability density of the YassoOT model parameters included in the package. 
Using this sample, it is possible to calculate the consequent uncertainty es- 
timates for the model predictions. Because of the availability of this sample, 
the uncertainty in the measurements, i.e. all the noise and variation not 
taken into account by the model, is inferenced directly into the predictions 

160 of the Yasso07-UI. 



4.1. Input and output information 

The Yasso07-UI works with input information that is easily accessible 
in practice. This information consists of: 1) the initial state of the soil in 
terms of the amount of ash free organic carbon in the soil and its chemical 
quality in terms of AWEN fractions; 2) the climatic conditions as monthly or 
annual precipitation and mean monthly or annual temperature; and 3) the 
estimated ash free organic carbon input into the soil in non- woody and woody 
litter, including the estimated size distribution of the woody litter - again 
in terms of the AWEN chemical fractions. The sequential extraction proce- 
dures to estimate carbon in the AWEN fractions are described in detail in 



e.g. iMcClaughertv et al.l (119851 ): iBerg et al.l fll991aU bl. ll993r ): lTrofvmow et al. 



(1l998l ). A list of chemical compositions of selected litters is provided at the 
YassoOT project website (Yasso07-UI manual). 



""This claim refers to a manuscript in preparation. However, some of these results are 
accessible through the YassoO? web site: www.environment.fi/syke/yasso. 
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The lack of accurate information on the initial state of the simulated 

175 soils is not limiting because given some constant litter input to the soil, the 
software can be used calculate the long-term stationary state of the soil to 
be used as the initial state. For modelling purposes, this state is completely 
defined using the total mass of OM in the soil and its division into the AWEN 
and H compartments in the model and by making the assumption that the 

180 litter input to the soils in terms of amount of carbon is equal to the output. 
The litter input to the soil is given to the software by defining the mag- 
nitude of the carbon flow into the soil in monthly of annual timescales. It 
can be given as a constant input or defined as a time-series describing e.g. 
the annual variations in the input or some long-term trend. The climatic 

186 conditions can also be given as some constant values or as time-series with 
monthly or yearly steps. This enables the user to investigate e.g. the effects 
climate change or short-scale changes in weather have on the soil carbon 
stocks and heterotrophic respiration. 

The Yasso07-UI presents the modelled estimates for the user in a variety 

190 of ways. The results are always available as raw numbers that describe the 
maximum a posteriori (MAP) estimate, the mean, standard deviation, and 
95% confidence intervals, or simply as a sample describing the probability 
density of the predictions. In addition, the software plots figures showing the 
amount of OM in each model compartment, the amount of OM in woody 

195 and non-woody litters, the total amount of OM, and the estimated respired 
CO2 function of time. 

4-2. Uncertainty in the Yasso07-UI 

The uncertainty estimates for model predictions are calculated in two 
phases. First, a sample is drawn from the joint probability density of the 

200 model parameters. When this sample is set large enough (preferrably at least 
few hundred), it represents the uncertainty in the parameter values of the 
model caused by the uncertainty and all the unmodelled excess noise in the 
measurements. In such cases, all the information in the measurements used 
to construct the model is inferenced to the model predictions as well. The 

205 sample size = Icorresponds to using only the MAP parameter estimates 
(Table [1]) without any uncertainty estimation. 

As a next step, the uncertainties in the initial state of the system (i.e. 
mass of OM and its chemical composition) and in the litter input are taken 
into account. The initial mass, and the percentage of mass in each compart- 

210 ment, are assumed to be Gaussian random variables. The means and the 
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standard deviations of these values are set by the user and the Yasso07-UI 
draws random values from these distributions. The litter input is treated 
similarly. 

The uncertainty estimation then proceeds as follows: 

1. A sample of size is drawn from the joint posterior density of the 
model parameters. 

2. For each parameter value 9i,i = 1, N, a random initial state is drawn 
from the PDF's of the initial state given by the user. 

3. For each parameter value 6i, a random input is drawn from the user 
defined PDF's of the input for the next time-step. 

4. The state of the system is calculated after the time-step for every 6i 
and corresponding initial state and input. 

5. This procedure is repeated by using the new state as the initial state 
of the next time-step. 

6. The basic estimates describing the probability density of the model 
predictions; such as mean, standard deviation, and 95% confidence 
intervals; are calculated for each time step in the output file. 

If the uncertainty of the input information is not known, we recommend 
the following settings. The uncertainty in the initial mass, namely its stan- 
dard deviation, should always be set to at least 10% based on the variations 
in the available data. This choise would account for the observed variability 
in the mass remai r iing measurements of approximately 10 - 15% of the initial 



mass fiBerg et all . Il993l : iGholz et all . l2000l : iTuomi et all . 120091 ). The uncer 



tainty in the initial chemical composition was found to be approximately 5% 
for each chemical component. This is approximately the amount of variation 



in the chemical composition of litter from a sinle plant species (iBerg et al. 
1991al lbh. 



For woody litter, the uncertainties set by the user should always be at 
least equal to th o se of n on- woody litter . Acco rding to the measurements of 
Palviainen et al.l (120041 ) ; IVavfova et al.l (120091 ) , the chemical composition of 
woody litter, in terms of AWEN compounds, does not vary more than that 
of non-woody litter. Also, the uncertainty in the litter input or initial state 
mass is not likely to be more than 10% for woody litter with diameter less 
than 10 cm. However, for woody litter with diameter greater than this, the 
uncertainty in the mass should be set to at least 15-20%. This is roughly 
the amount of variation in the measurements used to construct the woody 
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litter extension of the YassoOT model (ITuomi et al.l . l2010a[ ). Further, if e.g. 



the exact size distribution of the woody litter is not known, additional un- 
certainty should be added to the input and initial state mass estimates to 
account for all the actual uncertainty in the modelled system. 



5. Conclusions and discussion 

We have introduced a graphical user interface for the litter decomposi- 
tion and soil carbon model YassoOT. This user interface, Yasso07-UI, can 
be used to calculate point- and uncertainty estimates for litter decomposi- 
255 tion, soil carbon stocks, and CO2 emissions from soils in changing conditions 
with easily available input information. The user interface is applicable in a 
wide variety of climatic and environmental conditions because it takes into 
account all the information in the data sets used to build the YassoO? model 



(Tuomi et al.. 


2009 




2010a) 



260 The Yasso07-UI software can be used as a preliminary tool for assessing 
whether the YassoOT model can predict decomposition-related phenomena 
with respect to some specific set of measurements. These measurements could 
be e.g. litter mass-loss measurements, soil carbon stock measurements, het- 
erotrophic respiration measurements, or other measurements describing some 

265 features of the carbon cycle in soils. If the model appears to be consistent 
with the data, the YassoOT model can be easily implemented into different 
modelling purposes because its source code is freely available, the model has 
a simple structure, and the input information it requires is commonly readily 
available. 

270 Because of the nature of the measurements used to build the YassoOT 
model, there are cases where it is not known whether the model is applica- 
ble or not. These cases include swamps and marshlands, where the limited 
availability of free oxygen limits the decomposition heavily. This effect is 
not taken into account by our model. Also, the applicability of our model 

275 to modelling carbon stocks of mineral soils has not been tested throughly. 
Further, the model predictions regarding the decomposition of woody litter 
remain to be compared to wood decomposition data from temperate and 
tropical climates. We aim at developing the model further by extending its 
range of reliable applicability to these environmental and climatic conditions. 
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